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Chapter 9 TISSUE ENGINEERING 


9.1 Introduction 

Tissue Engineering (TE) is an interdisciplinary field dealing with the principles of engineering 
and fife sciences toward the development of biological substitutes that restore, maintain, or 
improve tissue function or a whole organ [52]. Currently, TE is emerging as an invaluable 
field of study and is one of the most promising fields for the next century. The ultimate goal 
is to develop powerful new therapies for biological substitutes that will successfully restore 
the structural and functional disorders of human and animal system. 

Available standard treatments for organ or tissue loss include transplantations of the 
same organ or organs from healthy donors, surgically reconstructed with synthetic im¬ 
plants in most of the cases, and the use of mechanical devices such as kidney dialysers, 
ventricular assist device (VAD) etc [14]. Of course, 
these treatments can save lives and decrease morbid¬ 
ity. However, each has its significant drawbacks and 
most are associated with substantially decreased qual¬ 
ity of patients fife. Although these procedures (treat¬ 
ments) temporarily save countless lives, unfortunately 
these techniques are not perfect. 

Tissue engineering procedures offer a new hope 
in modern medical science. By using patients own 
cells to produce bioengineered tissues and organs that 
function like their natural counterparts, TE may offer 
a better solution for the repair and reconstruction of 
damaged tissue and organs. At present, the only suc¬ 
cessfully engineered tissue that has been made avail¬ 
able to people is skin. Major success in this direction 
came in 1984 when Robert Langer, a chemical engi¬ 
neer at MIT and Joseph Vacanti, a medical surgeon 
at Harvard, first time jointly came up with an idea of 
duplication of bodys production tissues by developing 
a biodegradable scaffold on which skin cells could be 
grown using fibroblasts. After rigorous and extensive 
experiments with different techniques they discovered 
that collagen and fibroblasts can arrange themselves 
in a compact dense form under specific physical con¬ 
ditions that is exactly similar to the inner layer of the 
skin. Seeding keratinocytes on this layer leads to the formation of upper layer of the skin. 


Cartilage in 
external ear 


Cartilages in 
nose 


Cartilage 

In Kil«ow1«brd 

disc 


& 



Meniscus (padike 
cartilage in 
knee loinl) 




Putac 

symphysis 


Arbc Jar cartilage 
of a joint 


Key: 

B = Hyaline cartilages 
- Elastic cartilage* 
Fibrocartilages 


Fig 9.1: Possible locations of dif¬ 
ferent cartilages in human body 
( Source: [1]). 


1 



This is how the artificial skin was developed by TE methods. The important point of this ar¬ 
tificial skin is that it has no immunogenic cells, therefore question of rejection does not arise 
when implemented on patients body. Artificial skin are now made available commercially 
in large quantities. The next major step towards the development of artificial organ using 
tissue engineering approach is the cartilage. It is believed that the cartilage has relatively 
simple composition and structure. 

9.2 Cartilage Tissue Engineering 

9.2.1 Introduction 

Articular cartilage is a living tissue but of very 
limited capacity to repair itself. In fact, any dam¬ 
age to the articulating surface (due to injury, dis¬ 
ease, and/or genetic disorders) tends to spread 
allowing the bones to rub directly against each 
other and resulting in increasing joint pain and 
loss of joint movement. It can also lead to os¬ 
teoarthritis (a slow degeneration of articular car¬ 
tilage) (see Fig. 9.2). It is the most common form 
of arthritis in Australia, about 10% of the total 
population [61]. In United States alone, more 
than one million people annually suffer from os¬ 
teoarthritis, a degenerative joint disease caused 
by damage to articular cartilage [34], Damage to the articulating surface can occur due to 
compression of heavy load or when angular and/or shear forces are applied on it. With the 
prevalence of cartilage-related disease (e.g., Osteoarthritis), injuries and cost of treatment 
on the rise, the Osteoarthritis Research Society International (OARSI) urges new treatment 
options [34], An overview of a healthy and damaged knee is given in Fig. 9.2. 

The ability of cartilage to withstand enormous compressive loads (many times the body 
weight) without being crushed, is attributable to the multi-phasic nature of the tissue. Ma¬ 
ture articular cartilage contains approximately 5% of its volume as cells (also known as the 
growth cells or chondrocytes), and 95% as extracellular matrix (ECM). From an engineering 
standpoint, ECM is a porous, viscoelastic material consisting of three principal phases (a 
phase represents all of the chemical compositions with similar physical properties): 1) the 
solid phase constitutes about 30% of the total weight, and contains collagen (60%), proteo¬ 
glycans or PG (25%) and wide range of matrix proteins (15%), 2) the fluid phase (also known 
as the synovial fluid), which constitutes on average about 70% of the total weight, for normal 
tissues, and 3) an ion phase, consisting of many ionic species of dissolved electrolytes (e.g., 
Na + , Ca ++ , Cl - ). These constitute less than 1% of the net weight. The biofluid mechanics 
of cartilage growth and degradation relies on this multi-phasic nature of the tissue. 

9.2.2 Collagen 
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The collagen fibrils are densely packed polymer 
strands with a high resistance against fluid flow. 

They are the dominant structural components 
of the solid matrix and helps in retaining the 
shape of the cartilage when compressed. Col¬ 
lagens are basically made off three polypeptide 
chains arranged as triple helix long its length. 

There are two different types of collagens, known 
as fibrillar and non-fibrillar collagens. Depend¬ 
ing upon structural and functional characteris¬ 
tics collagen has been classified as 12 different 
categories. Hyaline cartilage (articular cartilage) 
mainly consists of type-II and type-XI fibrillar 
collagens where the percentage of type-II colla¬ 
gen is much higher. 

The architecture of articular cartilages connective tissues is very complicated. Collagen 
fibrils in the superficial-most region (known as the tangential zone) of the cartilage layer are 
densely packed and oriented parallel to the articular surface. This collagenous membrane 
has a relatively low PG content and a lower permeability to fluid flow, which is important in 
providing for a barrier of high resistance against fluid flow when cartilage is compressed. In 
the middle or transitional zone of the cartilage, the collagen fibers are larger and are either 
randomly or radially orientated. In the deepest zones, i.e., the zone nearest to the calcified 
cartilage and sub-chondral bone interface the collagen fibers are larger and form bundles 
that are oriented perpendicular to the calcified/bony interface. The way how the collagen 
and proteoglycan/glycosaminoglycans (GAG) molecules are arranged inside the articular 
cartilage, is shown in Fig. 9.3. 

9.2.3 Proteoglycan 

PG are macromolecules consisting of a protein core to 
which are attached short chains of negatively charged 
glycosaminoglycans (GAGs). The primary GAGs as¬ 
sociated with PGs in cartilage are chondroitin 4- 
sulfate, chondroitin 6-sulfate and keratan sulfate (Fig. 

9.4). It is available in all connective tissues but its 
structure varies with the type of tissue. Hyaline car¬ 
tilage contains mainly aggregating proteoglycans. 

Each PG-associated negative charge on the poly¬ 
mer chain (known as the fixed charge, FC), requires 
a mobile counter-ion (e.g., Na + ) dissolved within the 
interstitial fluid in the cartilage gel to maintain its 
electro-neutrality. This effect gives rise to an imbal¬ 
ance of mobile ions across the gel interface. The ex¬ 
cess of mobile ions colligatively yields a swelling pres¬ 
sure, known as the osmotic pressure, while the swelling pressure associated with FCs is 
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known as the Donnan pressure. Changes in this internal swelling pressure, arising from 
altered ion concentrations of the external bath, or changes in the fixed charges result in 
changes in tissue dimensions and hydration. However, this swelling pressure is balanced by 
tensile forces generated in the collagen network. This effect is due to the presence of covalent 
cross-links within the polymer matrix. An increase in cartilage tissue hydration, governed 
by the density and the nature of FCs on the PGs as well as the density of the mobile counter 
ions in the interstitial fluid are the earliest signs of articular cartilage degeneration during 
osteoarthritis. 

9.2.4 Chondrocyte 

Chondrocytes are mature cells found in cartilage. 

They make up the cellular matrix of cartilage, per¬ 
forming a number of functions within the tissue, in¬ 
cluding production and maintenance of the cartilagi¬ 
nous matrix and facilitating the exchange of fluids 
through the gelatinous layers. Fully mature chondro¬ 
cytes tend to be round, and they may cluster together 
in small groups within the network of the cartilage 
(Fig. 9.5). 

The progenitors of these cells arise in the bone 
marrow, in a form of stem cell. Stem cells are capa¬ 
ble of differentiating into several different types of cell, 
depending on the need. When they differentiate into 
cartilage cells, they start out as chondroblasts, ac¬ 
tively producing secretions of chondrin, the primary 
substance in cartilage, to build and repair the tissue. 

Once a chondroblast becomes totally surrounded, it 
is a mature chondrocyte. These cells can be found in Fig. 9.5: Light micrograph of hya- 
small gaps within the cartilage known as lacunae. line cartilage showing its chondro- 

Chondrocytes are not capable of cell division, cytes and organelles, lacunae and 
They can produce secretions to support and repair matrix ( Source: Wikipedia). 
the cartilage matrix, and as discussed above, the fa¬ 
cilitate the exchange of materials between the cartilage and the surrounding material. De¬ 
pending on what type of cartilage a cell is in, it may have a slightly different composition; 
elastic, hyaline, and fibrocartilage are all unique, designed to meet various needs within the 
body. 

9.3 Mathematical model of the ECM 

Mathematical models are essential for understanding systems (like cartilage) with many in¬ 
teracting components (e.g. the collagen network, proteoglycans, interstitial fluid and ionic 
solutions present in cartilage), since the complexity of the interactions taking place over¬ 
whelms the ability of the human intuition to make accurate predictions of their behavior. 
Comparisons of the models’ quantitative predictions with experimental results allows us 



4 



either to confirm that the biological hypotheses regarding these interactions (which are en¬ 
coded in the model) are correct or, equally usefully, identify that some of them are in error. 
Once thus validated, the model will ultimately help design treatments related to cartilage 
loss and osteoarthritis by predicting this outcome of proposed interventions, and focusing 
efforts on those which appear most promising. 

In this chapter, an account of the equilibrium for electro-chemical properties of poly(ethylene 
gylcol) gel (PEG) containing negatively charged chondrotin sulfate (ChS) strands, immersed 
in water with several different ions: H + , Na + , Cl - is demonstrated. This material closely 
mimics the electro-chemical properties of articular cartilage. We view this mixture as a 
multi-component material, including solvent particles, polymer and particles of several ionic 
species. The polymer is assumed to be made up of two types of monomers the uncharged 
units (PEG segment of the polymer) and those that are the charged ones (ChS segment of 
the polymer). Only the charged units participate in the binding reactions, carry a double 
negative charge and are denoted as M 2- . The positively charged ions in the solvent are 
Hydrogen (H + ) and Sodium (Na + ). The negatively charged ions could include Hydronium 
(OH - ) and Chloride (Cl - ). Because the negatively charged ions are assumed to be not in¬ 
volved in any binding reactions with the gel, acting only as counterions to positive charges, 
we identify these ions by the name Chloride. The binding reactions of the positively charged 
ions with the monomers are: 

k/7 k k/7 o 

(a) M 2- + H+ MH - , (6) M 2- + Na + MNa - , (c) MH - + H + MH 2 , 

k ~h k_ n k_^ 2 

(d) MNa - + Na + MNa 2 , (e) MH - + Na+ MHNa. (1) 

k—n 2 k —hn 

We assume that all the binding sites/charge sites are identical and the binding affinities for 
the different ions are different. The species M 2- , MH - , MNa - , MH 2 , MNa 2 and MHNa are 
different monomer species of the charged type, all of which move with the same polymer 
velocity. The ion species are freely diffusible, but because they are ions, their movement is 
restricted by the requirement to maintain electroneutrality. 

Suppose we have some volume V of a mixture comprised of k types of particles each 
with particle density (number of particles per unit volume) rq, and particle volumes zq, 
j = 1, • • • , k. From now on we will denote the quantities with subscript “1” related to PEG 
monomer species, subscript “2” related to ChS monomer species and subscript “3” related 
to solvent molecule. The subscript (p,s) denotes polymer and solvent phase, respectively. 
The volume fraction for each of the polymer species ($i,d 2 ) and the solvent (0 : > = 9 s ) are 
Oi = Viiii (i = 1,2,3) respectively and 9 P = 9\ + # 2 . Suppose the ChS monomers constitute a 
fraction ‘a’ of the total number of monomers, n p , (i.e., n 2 = an p , rq = (1 — a)n p ). Assuming 
that the ionic species do not contribute significantly to the volume (i.e., we take zq = 0, 
j = 4, • • • ,k), conservation of total volume implies 9\ + 0 2 + 0 3 = 1. Further, in subsequent 
calculations we make an assumption that the particle density of the ions is insignificantly 
small compared to the particle density of polymer and solvent, i.e. 4 rq <C n s + n p . 
Hence, 6 S = -f , 6 V = n , p , ob = — for j > 4 (assuming that the ions are dissolved 

in the solvent), are the polymer, solvent, and ion species per total solvent particle fractions, 
respectively. 
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9.3.1 Chemical potential 


This section outlines the calculation of the free energy and the chemical potential for this 
multi-species material. As above, we suppose there are k different kinds of particles, with 
particle volumes //, and particle numbers n,, i = 1,..., k. To calculate the chemical poten¬ 
tials, /ij, we use the Gibb’s free energy 

G = —k B TS + U + PV, (2) 

where U is internal energy, S is entropy, T is temperature, k B is Boltzmann’s constant, P 
is pressure, and V = ]>P ; UjUj is the total volume of the system. However, we assume that 
the volume occupied by species j = 4, • • • , k (which are the ions dissolved in the solvent) 
is small compared to that of the monomer plus solvent volumes v p n p + v s n s and so take 
V = v p n p + u s n s . The chemical potential are given by 


= —k B T - 


+ VjP — Ab T /T T zb-P, 


where //)’, /// are the contribution due to entropy and internal energy respectively. In the 
next two sections these two contributions to the chemical potential are outlined. 


9.3.2 Entropic Contributions to Chemical Potentials (//|) 

The entropy of the system is defined as 

S ^ ) TljUJi, (4) 

where ay is the entropy per particle for the i th particle. Using standard counting arguments 
given by Flory [23], for single-molecule species, 

ujj = - In (cf)j) j > 3 (5) 

The PEG and ChS chains exhibit permanent cross-link bonds (i.e., covalent bonds), and 
we calculate the per-particle entropy of these two species using the Doi’s rubber elasticity 
theory [20] 

=-^[(0i)' 2/3 - i], * = 1,2, (6) 

where i = 1,2 corresponds to PEG and ChS polymer species respectively. The corresponding 
particle fractions are (pi = (1 — a) (fi p and 4>2 = oU p . P and N t are the fractions representing 
the number of cross-linked monomers per total number of monomers in one chain, and 
the number of monomers per chain, respectively. Finally, the expressions for the entropic 
chemical potentials are given in the article by Sircar et al. [82] and are not included here for 
brevity. 
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9.3.3 Internal Energy Contribution to Chemical Potentials (pj) 

The internal energy consists of two contributions, long range electrostatic interactions and 
short range (nearest neighbor) interactions. The long range electrostatic interactions have 
energy 

U e = 'y ) ZjUj^e, (7) 

j 

where is the charge on the i th ionic species (z\ = z 2 = z p is the average charge per 
monomer), and $ e is the electric potential. 

To calculate the short range interaction energy for the polymer and solvent, we assume 
that each of the riT (= n p + n s ) particles have 2 neighboring interaction sites (called the 
coordination number). Of the total of n p monomers, k\n\ and fc 2 n 2 of them are the cross- 
linked PEG and ChS monomers, respectively. Since the cross-linked particles are pair-wise 
connected, we treat each of these pair as a single species. 

The different species with their pairwise interactions (i.e., the cross-linked species, x, 
the uncross-linked species, u, and the solvent particles, s) are shown in Fig. 9.6. Follow¬ 
ing assumptions are made while calculating the per-particle interaction energy, U: (1) all 
cross-links are assumed identical, (2) the cross-links are covalent permanent bonds, (3) the 
polymers have long chains and hence the end-effects are neglected, and (4) different polymer 
species have the same interaction energies if they are of the same type (e.g., uncross-linked 
PEG and uncross-linked ChS monomers have the same interaction energy). The cross-linked 
particle pair and the uncross-linked monomer have 2 . 2-6 and z-2 free interaction sites based 
on their position in the middle of the polymer chain, respectively (Fig. 9.6(a) and Fig. 
9.6(b)). The solvent particles have 2 free interaction sites (Fig. 9.6(c)). 
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Fig. 9.6: Pairwise interactions of the different species and their associated interaction ener¬ 
gies for: (a) cross-linked monomers, x, (b) uncross-linked monomers, u, and (c) solvent, s. 
Double solid lines denote cross-linking, single solid lines denote interactions within a polymer 
chain while the other types of interactions are denoted by dashed lines. The nearest neighbor 
interactions are assumed identical for both the PEG and ChS monomer species. 


The interaction energies for crosslinked, uncrosslinked polymer particles, and solvent 
particles are calculated using standard mean-held arguments. For example, the interaction 
energies for cross-linked particles, F x (Fig. 9.6a), are calculated by noting that there are 
either y- pairs of only PEG monomers, rj f pairs of only ChS monomers or ni ^ n ' 2 pairs of 
both PEG and ChS monomers, respectively. The corresponding probabilities of forming 
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cross-links are kf, k\ and k\k 2 . For each of these pairs, two neighboring sites are occupied by 
the covalent cross-link, four neighboring interactions are the monomer-monomer interactions 
within the polymer strands, and the remaining z — 6 interactions are distributed to monomer- 
monomer, with probability — and monomer-solvent, with probability —. Hence, 

^ + ~2^2 ■+-2-^1^2^ ( 2 e xx + ( 2 z — + 4 e w ), ( 8 ) 

where E us = e uu y + e us y. Similar arguments are made to calculate the interaction energies 
for uncrosslinked particles (F u ) and solvent particles (F s ). The total per particle interaction 
energy is 


U = 


1 


2(n m T n s ) 


{F x + F u + F s ) 


— ksTo 


X&ptps + I^O^S + ho ^P + o e w 


( 9 ) 


where the coefficients, X, ho)ho> are referred as the Flory interaction parameter and the 
chemical potentials of pure polymer and solvent species respectively [23]. These are functions 
of the particle fraction of the ChS strands (a), the interaction energy associated with the 
uncross-linked monomer-monomer interaction (ei), solvent-solvent interaction (£ 2 ), covalent 
cross-links interaction ( 63 ) and monomer-monomer bond within a polymer strand ( 64 ). The 
factor | in the per-particle interaction energy in Eqn. (|9j) is to correct for double counting. 
The corresponding contributions to chemical potentials (//j) are calculated via Eqn. <§• 


9.3.4 Interface conditions 

Under no-loading conditions, there is a free moving-edge to the gel on one side of which 
(inside the gel) 9 P = 6 ~, 6 S = 6 ~, and on the other side of which (outside the gel) 9 p = 0 , 
0 + = 1 , there are interface conditions, 

—Fp n = 0, (10) 

Z/p 

for the polymer, and 

—(/£-/O n = o, (11) 

Vs 

for the solvent, n is the normal to the free surface. The interface conditions are derived 
using the standard variational arguments to minimize the rate of work, which constitutes the 
viscous rate of energy dissipated within the polymer and the solvent, energy dissipation rate 
due to the drag between solvent and polymer and between solvent and ion species particles 
as well as the rate of work required against the chemical potential, py These interface 
conditions contain the contribution to the swelling pressure via the entropy, the Donnan 
potential and the net osmolarity (or the swelling pressure coming from the difference between 
the concentrations of ions dissolved in the gel and those dissolved in the bath, outside the 
gel). The derivation details can be seen in [46, 47, 82], 








9.3.5 Ionized-species chemistry 

We use the law of mass action to describe the chemical reaction. Under the assumption 
of fast chemistry and diffusion, the forward and the back reaction rates in Eqn. ([!]) are 
balanced. Those relations provide us with quasi-equilibrium concentration of the monomer 
species (e.g., m = [M 2- ],r/ = [MH“],r = [MNa - ],w = [MH 2 ],r = [MNa 2 ] and q = [MHNa]). 
Similarly, the quasi-equilibrium concentration of the ion species ( n = [Na + ], h = [H + ] and 
ci = [Cl - ]) reduces into 

C = C b e~ zc *% (12) 

where C = h,n, ci, z n = Zh = 1 and z C[ = — 1 and the subscript ‘b’ denotes the corresponding 
bath concentrations (or concentrations outside the gel). The electrostatic potential, T e , is 
determined by the electro-neutrality constraint inside the gel, namely, 


(n + h — ci)9 s + z p rriT = 0, 


(13) 


where z p is the average charge per monomer. This charge depends on the residual charge of 
the unbound and the bound ChS monomers, z p , given by 


z p rriT = —(2m + y + v). 


(14) 


From Eqns (13, 14), we can calculate the Donnan potential, z p ty e 
outside the gel (since there is no polymer in that region). 


This potential is zero 


9.4 Experimental methods 

In vitro experiments provide the most intuitive mechanism for cartilage swelling-deswelling 
kinetics. Recent experiments have confirmed that changes in external ion concentration 
and pH affect not only the mechanical properties of cartilage, but also the biosynthetic 
capabilities of chondrocytes. The experimental outcomes of Beekman et al. [7], Freed et al. 
[25], Kisiday et al. [49], Regan et al. [73], Stading and Langer [84], Vunjak-Novakivic et al. 
[96] and Wilson et al. [99], all share a common feature: collagen accumulation in engineered 
cartilage tend toward a constant value (steady state) with respect to time. A similar type 
of result was obtained by Buschmann [16] for GAG molecule accumulation in engineered 
cartilage constructs. 

According to the results on DNA content of engi¬ 
neered cartilage constructs obtained by Wilson et al. 

[99] (Fig. 9.7), there is an increasing trend of cell pop¬ 
ulation growth that means the cell population has no 
control (i.e., it keeps on increasing unboundedly), but 
this is not the case in nature. Hence, overall experi¬ 
mental observations suggests that ECM deposition in 
cell-polymer construct of engineered cartilage heading 
towards a steady state but there is no concrete evi¬ 
dence that chondrocytes population does the same. 

Our experimental validation of the model, de¬ 
scribed previously, is performed by running a subset 



Fig 9.7: DNA concentration Source: 
[99] 
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of in vitro experiments in tandem with the proposed 
computational experiments described below. We con¬ 
struct ionic gels by co-polymerizing methacrylated chondroitin sulfate and poly(ethylene gly¬ 
col) dimethacrylate in deionized water (distilled water). PEG was chosen since it can be used 
as a neutral, cross-linked, 3-D scaffold, a system that promotes the deposition of proteogly¬ 
cans and collagen molecules and emulates the mechanical strength, load bearing capabilities 
and resilience of cartilage tissue. The ChS component serves two purposes: to mimic the 
biochemical environment of cartilage (since it is the main component of proteoglycans) and 
to introduce fixed negative charges into the network. These polymer gels are then allowed 
to reach equilibrium in deionized water (e.g., 24 hours) and then placed into appropriate 
solution (e.g., with variable pH, salt concentration, etc.) to reach a new equilibrium. 

To quantify the local configuration of the polymer gels, the counter-ion concentrations and 
the osmolarity inside these gels, the following in-vitro experiments were performed and the 
results were utilized to validate the outcome of numerical simulations: (i) Vary the particle 
fraction of PEG:ChS in a salt-free, neutral bath (H*, = 10 -7 M), where a, the particle fraction 
of ChS, is: 0.1, 0.2, 0.3, 0.4, 0.5 and 1. This experiment is performed to determine the role 
of fixed charges. Since the ChS monomer volume and the density are constant, increasing 
the particle fraction of ChS directly corresponds to increasing its weight fraction, (ii) Vary 
the bath concentration of the monovalent salt within a range of concentrations, 0.1-1 M, and 
(iii) repeat experiments (i) and (ii) in a salt-free bath and variable pH for neutral hydrogels 
(i.e., a = 0) as well as for ionic gels (i.e., a > 0). The pH is varied in the range 3-8 to 
ensure that the model accurately captures the local variables (i.e., the volume fractions and 
the osmolarity), within the physiological range [57]. Further, the above mentioned three 
sets of experiments are repeated for gels with different cross-link fractions. Gels of varying 
cross-link density are achieved by changing the total monomer concentration in solution prior 
to polymerization (e.g., 10 wt %, 20 wt %, and 30 wt % gel). Polymer volume fraction is 
determined from equilibrium swollen volume and dry polymer volume measurements and 
compared to the numerically determined value. The initial volume of the dry polymer, V*, 
was determined from the initial mass, up, using a weighted density of PEG and ChS, as 


follows 


Vi = 0.1 


rrti ■ %PEG mi ■ %ChS 
Ppeg Pens 


(15) 


The mass of the gel at equilibrium was measured and used to determine the swollen volume, 
V swollen- It was assumed that all polymer was incorporated in the gel, which was confirmed 
for chondroitin sulfate, and therefore a change in mass was solely due to the a change in the 
water content. Hence, for polymers with a 10 wt % dry weight in solvent, 


/ rn,j • %PEG m j • %ChS \ mswollen 0.1 m, 

Vswollen = 0.1-1-1 H-(16) 

V Ppeg PChS 7 Psolvent 

The initial mass, m*, was the mass of the gel immediately following polymerization and before 
being placed in a solvent. The swollen mass, m swo u en , is the mass of the gel at equilibrium. 
The volume fraction of the polymer was calculated as the volume of the dry polymer (i.e. 
before swelling) divided by the volume of the swollen gel, Vi/V swo u en . The kinetic chains 
associated with methacrylate group (-MA) represents less than one percent of the total gel 
volume and therefore are not considered in the volume calculations. 
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PEG (i = 1) 

ChS (* = 2) 

Units 

Density (p/) 

1.07 

1.001 

g/mL 

Molecular weight (Mj) 

4600 

48700 

g/mol 

Repeat unit per chain (AT*) 

102 

86 

- 

Cross-link fraction (k/) 

0.25, 0.5, 0.75 

0.25, 0.35, 0.45 

- 

Hildebrand solubility (Si) 

17.39 

5.19 

MPa 1 / 2 

Hildebrand solubility for wa 

ter (<S W ) 

48.07 

MPa 1 / 2 


Table 9.1: Parameters common to all the numerical results. The reference temperature for 
the solubility parameters is fixed at T 0 = 298K. 


9.4.1 Comparison of model outcomes with experiments: parameter estimation 


The values of the parameters used in numerical calculations are listed in Table 9.1. The 
constants in the model are the monomer volumes, ui,u 2 , the coordination number of the 
PEG-ChS polymer lattice, z, and the nearest neighbor interaction energies, e* (§9.5.2). The 
undetermined constants are the binding affinities of the polymer gel with Sodium ( K n ) and 
hydrogen (K n ). 

The monomer volumes are found from the density and the molecular weight information 
(//, = Mi/{pi * Na))- The cross-linked PEG-ChS matrix has a 3-D configuration, which 
suggests that we choose the coordination number, z = 6, which mimicks a 3-D cubic lattice. 
The standard free energies, an d ksTop,® and the interaction energies, e, are found 

from the Hildebrand solubility data, Si [?]. The standard free energy is the energy of all the 
interactions between the molecule and its neighbors in a pure state that have to be disrupted 
to remove the molecule from the pure state. 

The molecular mass and the density of PEG were fixed at 4600 g/mol and 1.07 g/mL, 
while that of ChS were selected at 48700 g/mol and 1.001 g/mL, respectively. These values 
give the monomer volumes of the PEG chains as zq = 7.14 x 10 -21 cm 3 , and that of ChS 
chains as zq = 8.08 x 10 -20 cm 3 . The density of the solvent, p so ivent, was assumed to be 
1.00 g/ml and the volume of 1 molecule of water at reference temperature, To = 298K, is 
u w = 2 x 10 -23 cm 3 . The Hildebrand solubility parameters for pure species (values given 
in Table ??) and the monomer volumes are used to calculate the interaction energies, e t 
(i=l,...,4). Using the solubility relations [82], these values are fixed at e\ = 0.0, = 3.74 ,63 = 

9 . 58,64 = —56.17. The reference temperature of the experiments was fixed at T 0 = 298K. 
The undetermined parameters, namely the binding affinity of the polymer gel with various 
cations (Kh, K n ) are computed by minimizing a nonlinear least-square difference function 
between the experimentally observed values of the equilibrium volume fraction, 6 P , and 
the corresponding model output, implemented via the MATLAB non-linear least-square 
minimization function lsqnonlin. These values are found as logio(iP n ) = —2.17, logio(AT) = 
-3.65. 


Fig. 9.8 presents the closeness of fit between the model (highlighted by the solid lines) and 
the sample averaged equilibrium data-points at different salt concentrations in a neutral pH 
sample (Fig. 9.8a) as well as in a salt-free, variable pH solution (Fig 9.8b). The experimental 
values were noted at time t — 24 hrs. Each point represents the average of four samples. 
The values are noted in identical conditions with the upper and the lower limits in the error 
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Fig 9.8: Experimental data of equilibrium volume-fraction of the PEG-ChS gel vs. (a) 
different NaCl concentration in the bath (in mol/lt) and pH = 7.0, and (b) different pH and 
zero salt concentration. The gel sample is 10 wt % in solvent. The cross-link fraction of the 
PEG/ChS polymer pairs are k\ = 0.25 and k 2 = 0.25, respectively. Notice the closeness of 
fit between the model output (i.e., solid lines) and the experimental data points. 


bar representing maximum and the minimum variation from the average, respectively. 

9.5 Effect of structural and environmental fluctuations on equi¬ 
librium ECM configuration 

Articular cartilage is composed of highly negatively charged aggrecan aggegates dispersed 
in a crosslinked collagen network. This chemical make-up gives rise to an unusually high 
concentration of mobile-counter ions (e.g., Na+) in the interstitial fluid of the tissue. As 
a result, the cartilage tissue experiences a complex set of electromechanical cues (signals) 
arising from fixed negative charges (aggrecan), and mobile-counter ions, in equilibrium. Un¬ 
derstanding and recapitulating the electromechanical cues in a tissue engineering approach 
may be key to bioengineering functionally competent cartilage. 

9.5.1 Effects of changes in ionization, cross-links 

Fig. 9a,b,c present the equilibrium volume fraction, corresponding Donnan potential and 
the net-osmolarity versus the particle fraction of ChS, «, respectively. The various curves 
represent gels with different cross-link fractions, (ki,k 2 , Eqn. §• It can be shown that 
the average charge per monomer (or the average ionization) increases in proportion to the 
particle fraction, a. The gel is dissolved in salt-free, neutral water. Any far field boundary 
effects are neglected in our numerical simulations, by assuming infinite bath conditions. The 
negative charges on the gel causes water to dissociate into hydrogen (H + ) and hydronium 
(OH - ) ions. The free H + ions bind with the monomers. 

At low ionization (or low particle fraction values, a), the physically relevant solutions 
swell as a function of increasing a. This is due to the increased swelling effect via the 


12 








Fig. 9.9: (a) Equilibrium polymer volume-fraction, (b) Donnan potential, and (c) Net- 
Osmolarity in neutral salt-free, bath conditions vs. a, the particle fraction of chondrotin 
sulfate component for gel solutions with different cross-link fractions. 
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Donnan potential (e.g., look at the trend of the volume fraction and the Donnan potential 
in the range, 0 < a < 0.9). The physically relevant solutions are the largest and smallest 
ones, if there are three solutions, because the intermediate solution is unstable. Another 
trend is that the gel deswells versus increasing a, at high cross-link fraction and particle 
fraction values (the curve denoted by CL 3 in the range a > 0.9, figure). This observation is 
corroborated by experimental results that a highly cross-linked network provides a barrier 
against the Donnan potential and deswells at higher ionization levels [48,67]. 

9.5.2 Effects of changes in the bath salt concentration 

Figure 9.10 a,b,c highlight the corresponding equilibrium volume fraction, Donnan potential 
and the net-osmolarity versus the bath salt concentration, [NaCl]&, for polymer gels immersed 
in a neutral water solution, respectively. These gels have variable cross-link fractions but a 
fixed ChS particle fraction, a = 0.1. Usual features emerge at high concentrations ([NaCl]& > 
10 -3 M) and low concentrations ([NaCl]/, < 10 -8 M) of salt. At higher salt concentrations, 
the gel deswells and the deswelling is more prominent for specimen with greater cross-link 
fraction. This is readily understood since at higher salt concentrations, more dissolved Na + 
ions are furnished which can bind with the negatively charged gel. As a result, the primary 
ingredient of swelling pressure, i.e., the Donnan potential, is negligible (Fig. 9.10b). High 
cross-link density within the gel resists the swelling pressure which further enhances the 
deswelled state of the gel. At lower concentrations of salt, there any hardly any Na + (or H + ) 
ions readily available to bind with the gel. The Donnan potential is nearly constant and this 
feature leads to a fixed equilibrium state of the gel ( 9 P ~ 0.65). 

For intermediate salt concentrations (10 -8 < [NaCl]& < 10 -3 M), a new swelling-deswelling 
feature emerges, namely, there is a swelling, either gradual or through a phase transition, 
as the salt concentration in the bath increases. The explanation for these swelling feature 
is that there is a complicated interplay between ionization (which promotes swelling) and 
increases in bath concentration of salt which promote deswelling [82], 

9.5.3 Effects of changes in the bath pH 

Two sets of numerical experiments are performed to predict the equilibrium volume fraction 
of the charged gel. In the first experiment, cross-link fraction of the polymer matrix is 
fixed, ki — k -2 — 0.25, while the ChS particle fraction is varied. Overall, it is found that 
increasing the bath concentration of hydrogen leads to deswelling. The reason is that a 
decrease in pH leads to an increase in the dissolved cations which readily bind with the gel, 
thereby decreasing the gel ionization and the swelling effect via Donnan potential. In this 
case, swelling is predominantly driven by the ionization which is induced by the Donnan 
potential. The pressure created by the Donnan potential is about an order of magnitude 
greater than potential induced osmotic pressure (e.g., compare the order of magnitude of the 
curves in Fig. 9.11b versus Fig. 9.11c). 

Another set of numerical simulations is carried to determine the equilibrium solutions of 
gels having a fixed ChS particle fraction but a variable cross-link fraction. It is found that 
gels progressively preferred a deswelled state in the order of increasing cross-link fraction 
values, consistent with our observations in 9.5.2. In summary, for the values of parameters 


14 






Fig. 9.10: (a) Equilibrium polymer volume-fraction, (b) Donnan potential, and (c) Net- 
Osmolarity, vs. NaCl-concentration (in mol/lt.) in the bath. Different curves represent 
samples a fixed ChS particle fraction, a = 0.1, and variable cross-link fractions. Note the 
match in the equilibrium volume-fraction and the experimental data (Fig. 9.8a) in the range, 
—4 < log w [NaCl]b < 0. 
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examined, the gel swells predominantly due to effects of increasing Donnan potential and 
deswells due to the combination of increasing bath concentration of cations and increasing 
covalent cross-links within the polymer matrix. 





(k) 


Fig. 9.11: (a) Equilibrium polymer volume-fraction, (b) Donnan potential, and (c) Net- 
Osmolarity vs. hydrogen bath concentration H for a gel with fixed cross-link fraction 
ki = k 2 = 0.25 and variable composition of ChS. Note the match in the equilibrium volume- 
fraction and the experimental data (Fig. 9.8a) in the range, —8 < log 10 [i7](, < —3. 


9.6 Concluding remarks 

In this chapter a description of the makeup of the cartilage tissues is given. We then pro¬ 
vided a multi-phase, multi-species model to shed light on the swelling / de-swelling mech¬ 
anism for ionic gels with covalent cross-links; as well as provided a detailed description of 
the experimental methods to engineer these biogels. These polymeric biogels mimic the 
electro-chemical environment of articular cartilage. Finally, we have highlighted the effect 
of structural and environmental uctuations on the equilibrium tissue configuration. 
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Current research includes determining the counter-ion concentration and local osmolarity 
in dynamically loaded gels, determining the effect of the environmental fluctuations on the 
cartilage cell response and the tissue development. A synchronous development of biome¬ 
chanical models and numerical schemes are following. Readers are requested to look at the 
research work done by Stephanie Bryant et al. [13,95] on the experimental side and those 
by Van Mou et al. [51,65,66] on the numerical end. Ultimately, in the long term, the scien¬ 
tists and the tissue engineers expect to answer the fundamental question through research in 
this area, how important are electro-mechano-chemical cues in guiding cellular response and 
cartilage tissue growth? Information gleaned from models and experiments will help under¬ 
stand the pathway to cartilage degeneration, repair strategies, maintaining healthy joints via 
delayed tissue aging and developing superior mechano-chemically resistant biogels as joint 
replacements. 
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